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Abstract 

We study the matrix model for N M2-branes wrapping a Lens space L(p, 1) = 
S 3 /Z, p . Lhis arises from localization of the partition function of the ABJM 
theory, and has some novel features compared with the case of a three-sphere, 
including a sum over flat connections and a potential that depends non-trivially 
on p. We study the matrix model both numerically and analytically in the large 
N limit, finding that a certain family of p flat connections give an equal dominant 
contribution. At large N we find the same eigenvalue distribution for all p, and 
show that the free energy is simply 1 jp times the free energy on a three-sphere, 
in agreement with gravity dual expectations. 
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1 Introduction and summary 



One of the most exciting features of supersymmetric gauge theories is that one can 
compute certain protected quantities exactly. The most fundamental of these quantities 
is the partition function. In recent years localization techniques have been developed 
that allow the computation of partition functions for supersymmetric gauge theories 
in different dimensions, and on different backgrounds, starting with the work of [HI2]. 
These results have also led to new tests of conjectured dualities between theories. 

Three-dimensional M = 2 gauge theories form a particularly fertile ground in which 
to develop these ideas. The simplest compact manifold on which one can define d = 3 
supersymmetric field theories is the round three-sphere, originally studied in [2]. This 
generalizes to other three- manifolds M3, including S 1 x S 2 and certain one-parameter 
families of squashed three-spheres [31 H] . In this paper we study M = 2 theories on the 
Lens spaces L(p, 1), which are free quotients of S 3 by Z p . In the first part of the paper 
we derive a formula for the full localized partition function of a three-dimensional M = 2 
U(N) Chern-Simons-matter theory on such a Lens space. Several of the ingredients 
have already appeared in previous papers, including [SI El El- The partition function 
reduces to a matrix model integral, where the potential function depends non-trivially 
on p. 

In the second part of the paper we consider the partition function in the large N limit, 
keeping the Chern-Simons levels fixed. The motivation for this is that, for appropriate 
matter content, one expects to be able to reproduce these results from a dual M-theory 
gravity computation. In particular, we focus on the low energy effective theory on N 
M2-branes, described by the AB JM theory [8] . A new feature that arises when M3 has 
non-trivial fundamental group is that one must sum over different topological sectors in 
the partition function. In the present case, different sectors are labelled by a diagonal 
N x N matrix with entries in Z p . In the large N limit, we show that one can in fact 
focus on the contribution from matrices proportional to the identity. This drastically 
simplifies the analysis. 

At large N we may use a saddle point approximation to the matrix model, follow- 
ing [9]. The leading contribution to the free energy arises from a specific eigenvalue 
distribution. In order to gain some intuition we study this distribution numerically, 
for a number of values of p. These numerical results lead to a simple ansatz for the 
eigenvalue distribution at large N. We then use this ansatz to obtain analytic results 
for the free energy, as well as for the eigenvalue distribution and corresponding density. 
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We find that the eigenvalue behaviour is in fact independent of p, with a free energy 
that is simply 1/p times the free energy on a three-sphere. This is in agreement with 
gravity dual expectations, where L(p, 1) arises as the conformal boundary of AdS^Zp. 

The organization of this paper is as follows. In section [2] we derive the full localized 
partition function for a three-dimensional M = 2 Chern-Simons-matter theory on the 
Lens space L(p, 1). Section [3] contains the numerical results, which are compared to 
corresponding analytic results in section |4| We mention some open problems in the 
outlook section [5j Finally, some technical results are relegated to the appendices. 

2 The localized partition function on L(p, 1 



In this section we derive a formula (2.25) for the full localized partition function of 



a three-dimensional M = 2 Chern-Simons-matter theory on the Lens space L(p, 1) = 
S 3 /7* p . We then specialize to the ABJM theory [5] on N M2-branes of interest. 

2.1 The Lens space L(p, 1) 

The Lens space L(p, 1) = S 3 /7* p is a certain freely-acting quotient of the round S 3 
by a group of order p. Regarding S 3 as a unit sphere in Euclidean C 2 , with complex 
coordinates zx,z 2 , the Z p action is generated by 

C 2 3 (z 1 ,z 2 ) H> (u p z 1 ,u~ 1 z 2 ) , (2.1) 

where 

Up = e 2ni / p is a primitive p-th root of unity. Notice this is simply a Z p quotient 
along the S 1 fibre of the Hopf fibration: S 1 ^ S 3 — > S 2 . 



Quotienting S 3 by the free action (2.1 ) leads to a smooth three-manifold L(p, 1) with 
iTi(L(p, 1)) = Z p . There are then p topologically inequivalent complex line bundles L 
over L(p, 1), labelled by their first Chern class c\(L) £ H 2 (L(p,l);Zj) = 7L V . Each L 
admits a flat U(l) connection, which plays an important role in studying gauge theory 
on L(p,l). For example, rather than complex- valued functions on L(p, 1), it will be 
important to consider more generally sections of L. 

Concretely, one can construct such sections as certain projections of functions on 
the covering space S 3 . For example, we may expand complex-valued functions on S 3 
in terms of hyperspherical harmonics 

Y^ n {e,M = », m ,n(^)e im ^, (2.2) 
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where (6, <f>, if>) are standard Euler angles on S 3 . Here i G Z>o while m,n G {- 



2 ' 



+ 



2 



}, with £ labelling the 



2' 2. 



spin representation of SU(2)l X SU{2)r acting on 



S* 3 . For example, Ye )mjn has eigenvalue — + 2) under the Laplacian. 

If Zi = Tie 1LPl denote polar coordinates on each copy of C in C 2 , then <f> = <pi + <p 2 , 
if) = <fi — (f2, and in terms of Euler angles the generator (2.1) thus acts as 

Air 



if) i — y if) -\- 



P 



(2.3) 



The complex-valued functions on L(p, 1) are precisely the Z p -invariant functions on 
S 3 , and thus (2.3) and (2.2) imply that functions on L(p, 1) are spanned by the modes 

2n = mod p . (2.4) 



(2.2) satisfying 



More generally, sections of L are spanned by the modes (2.2) satisfying 

c\(L) mod p , 



2n 



(2.5) 

This follows since the holon- 



where we are using the isomorphism H 2 (L(p, 1); Z) = Z p 
omy of the flat connection A on L around the generator 7 of 7r 1 (L(p, 1)) is 



exp 



A 



^vrici (L)/p 



(2.6) 



Here 7 is represented by a circle fibre in S 1 ■=->■ L(p, 1) — )■ S* 2 . Sections of L must then 
also pick up this phase around 7. 

Another issue, important for considering supersymmetric field theories, concerns the 
Killing spinors. The 4 Killing spinors on S 3 transform in the (2,1) © (1,2) repre- 
sentation of SU{2)l x SU{2)r = Spin(4). The spinor used for localization in [2] is 
in the (2,1) representation. In this language the Z p action in (2.3) is contained in 
U(1)r C SU{2) l x SU(2) R , and hence the Killing spinor used for localization on S 3 
projects down to a Killing spinor on L(p, 1). 

We pause here to make some comments on more general Lens spaces L(p,q), with 
q > 1. These are defined as the free quotient of S 3 C C 2 by the action 

C 2 3 ( Zl , Z2 ) _> (rfzuu; 1 ^) , (2.7) 



with q relatively prime to p. The Z p action (2.3) now becomes 

1) 



1 ^ 



+ 



2n(q 



if) 1 — y if) -\- 



P 

2vr(g + l) 



P 



(2.8) 



on the Euler angles. In particular, there is therefore no invariant spinor, unless q — 1 
or q = p — 1. This makes the treatment for this case more involved, and we therefore 
leave it for future work. 



2.2 Localization of the path integral 

The localization of the path integral on L(p, 1) is very similar to the original computa- 
tion for S 3 in [2\. Indeed, locally the spaces are identical, so one just needs to keep track 
of how global differences affect formulae. For example, for a U (N) gauge theory the 
path integral still localizes onto flat connections A, but on L(p, 1) there are non-trivial 
flat connections that one must then sum over. A flat U(N) connection on a manifold 
M is determined by its holonomies, which define a homomorphism g : 7i"i(M) — > U(N). 
Gauge transformations act by conjugation, so that flat U(N) connections are in 1-1 
correspondence with 

Hom(vri(M) U(N)) /conjugation . (2.9) 

Since 7Ti(L(p,l)) = Z p , specifying a flat connection is equivalent to specifying the 
holonomy around the generator 7 of 7Ti(L(p, 1)) = Z p 

hol 7 (A) = diag«\...,<») , (2.10) 

where < m, < p, and % = 1, . . . , N runs over the generators of the Cartan U(1) N 
subgroup of U(N). Here we order mi < 1712 < ... < rriN (conjugation permutes 
the entries). The localized path integral will then give a sum over topological sectors 
m = diag(mi, . . . , m N ), for each U(N) gauge group. 

Apart from this, as for S 3 all fields localize to zero except for the D-term and scalar 
a in the M = 2 vector multiplet, which are related via 

D — —a . (2.11) 

The scalar a must be covariant constant. Writing the flat gauge field defined by m 
as A m = —ig^dgm, this implies that a = g^^oOm where a is a constant N x N 
Hermitian matrix satisfying 

K,m] = 0. (2.12) 

For a Chern-Simons gauge theory, this saddle point solution gives a standard classical 
contribution to the saddle point approximation of the path integral 

i7rfc 



exp [-Sdassicai(o"o)] = exp 



p _Tr (°o) 



(2.13) 
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coming from the supersymmetric completion of the Chern-Simons interaction, evalu- 



ated on (2.11). Here k G Z is the Chern-Simons level. The p-dependence in (2.13) 



simply arises because Vo\(L(p, 1)) = Vo\(S 3 )/p. The path integral then reduces to a 
matrix integral over o"o, as well as the discrete sum over m labelling flat U(N) gauge 
fields. One must also include the Chern-Simons action for the flat gauge field: 



ex P [-Scs(A)] = exp 



exp 



ik 

Atx 
mk 
P 



Tr 



Tr (m 2 "l 



A A dA + -A 3 
3 



(2.14) 



One computes (2.14) in a standard way: choose a four-manifold M4 with boundary 
9M4 = L(p, 1), and an extension of the bundle and (flat) connection A on L(p, 1) 
to corresponding data over M4. The Chern-Simons action is then in fact defined as 
— (ifc/47r) J M Tr (F A F), which can be shown to be independent of choices, modulo 
27ri. For example, in the present case one can take M 4 = total space of 0{p) — > CP 1 , 
and note that the restriction map Z = if 2 (M 4 ;Z) — > H 2 (L(p,l)]Z*) = Z p is simply 
reduction mod p. 

Having summarized the localization, we next turn to the effect on the one-loop 
contributions around the saddle points specified by (<7o,m). Due to the remarks in 



section 24^ the spectra of operators that contribute to the one- loop determinants reduce 
to an appropriate projection of the full spectra on S 3 . 

2.3 Matter multiplet 

We consider here the contribution of a chiral matter field in the representation 1Z of 
the gauge group, to the one-loop determinant around the classical background labelled 
by (00, m). We denote the R-charge of $ as A = A($) - the canonical value is A = | 
- and the weights of 1Z by p. 

The bosonic contribution to the one-loop determinant is then [21 [101 [H] 

det, /2 , m (D boson ) = £(£ + 2) - 4m(l - A) + A 2 + p(a ) 2 , (2.15) 

where £, m label the same quantum numbers as in section 2.1[ so that i > and 



m G {- 



2 ; 



+ 1, 



• 2 



}. In particular, the l[l + 2) term simply comes from the 



eigenvalue under (minus) the scalar Laplacian. On S 3 there are I + 1 such modes, 
labelled by the quantum number n E { — f , — 1 + 1,...,|}, while on L(p, 1) we should 
keep only those modes satisfying 

2n = p(m) mod p , (2-16) 
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as follows from (2.5) 



The fermionic contribution to the one-loop determinant is also given by (2.15), 
but now with m G { — |, — | + 1, . . . , | — 1}, and with an additional contribution of 
{—l) e {£ + A + ip(o"o)) (£ + 2 — A + ip(a )). Again, on S 3 there is a degeneracy of £ + 1, 



labelled by n, while on L(p, 1) we should keep only those modes satisfying (2.16). 



Since the one-loop determinant is a ratio of fermionic and bosonic determinants, we 



thus see that for fixed £,m and m, for every choice of n satisfying (2.16) the contri- 



butions from fermionic and bosonic determinants will cancel, except for the "missing" 
fermionic mode with m — | - this remains uncancelled in the bosonic determinant. 
We thus conclude that, for fixed £, we have 

fermion 

) , (£ + A + ip{a )) {£ + 2 - A + ip{a )) 



boson 

) v ' £(£ + 2) - 2£(1 - A) + A 2 + p(a ) 



_x 2 

+ 2 - A + ip(o-o) 



(2.17) 



v ' £ + A-ip(a ) ' 
where the degeneracy is the number of half-integers n E { — |, — 1 + 1,...,|} satisfying 

2n = p(m) mod p . (2-18) 

This degeneracy was denoted by N p (£) in [?]. Thus in total 



2.4 Vector multiplet 

The analysis for the one-loop contribution of the vector multiplet is very similar. Here 
it is more convenient to follow the analysis in [Rj, where rather than working out the 
full spectrum, most of which then cancels in the ratio of determinants, instead one 
isolates the uncancelled modes from the outset. These are precisely the eigenmodes 
which are not paired with a superpartner. We shall refer to these as the uncancelled 
modes. 

The uncancelled gaugino modes on S 3 have eigenvalues 

p = n 1 + n 2 + ia(a ) (2.20) 

under the relevant Dirac operator, where rij denote the charges under d 9i , where recall 
that ipi, f 2 are azimuthal angles on each copy of C in C 2 D S 3 . The normalizable 
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modes are {ni,n 2 > 0} \ {(rii,n 2 ) = (0,0)}. The corresponding charges under d^, 



are then |(m — n 2 ), + n 2 ), respectively (see just before equation (2.3)), so that 
the projection condition becomes 



n 2 + a(m) mod p . 



(2.21) 



The uncancelled transverse vector modes also have eigenvalues of the form (2.20) 
except now ni,n 2 < —1. Thus 



^rsSpKm) = n 



n 



n l i n 2 ^ 

(ni,n 2 ) jt (0,0) 
?q ri2 + a(m) mod p 



[ni +n 2 + ia(cro)] 



n 



ni = ri2 — a(m) mod p 



[-ni - n 2 + ia(er )] 



(2.22) 



We then rewrite this as a product over only the positive roots a > 0, while at the same 
time multiplying by the same expression with a — > —a. In doing this, one sees that 
all the terms in the numerators and denominators cancel, except for the numerator 
contributions of {n\ = 0,n 2 > 1}, {ni > l,n 2 = 0}, which are left uncancelled. We 
are thus left with 



^r C io O o r p(^0,m) 



n r , v . , ■ / m rir>o[«( m ) +pr-ia(a )] 



n 

o>0 



n r>0 [a(m) +pr +ia(a )} 



r>l 



{ia(a )Y 



»(m),0 



Y[ [-Qt(m) + pr - ia(a )) 



r>l 



. (2.23) 



Notice that here, in a slight abuse of notation, we have assumed that < ct(m) < p. 
The last equation may then be rewritten 



n 

a>0 



r=l 



[aim) — ia{<jQ 



»n 



r=l 



(a(cr ) + ia(m)Y 



X 



(a(m) + ia(a )) 

r=l 



1 + 



(a(g ) - ia(m)) z 



Zeta function regularizing n^Li(P r ) 4 ^ (2tt) 2 /p 2 and using the infinite product for- 
mula for sinh, we obtain [51 U\ 



= n 



2 sinh 



a>0 



£ (a(cr ) + ia(m)) 
(a(cr 



2 sinh 



I (o(^o) - ia(m)) 



. (2.24) 
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2.5 Partition function 



Putting everything together from the previous sections, we arrive at the final formula 
for the localized partition function of an M = 2 Chern-Simons-matter theory on the 
Lens space L{p, 1) 

mk 
. P 



Z = J danexp 



(Tr(a 2 )-Tr(m 2 )) 



^vector / 
-^l-loopl 



a ,m)Z 1 m _T o t o e pKm),(2.25) 



where the one-loop vector and matter contributions are given by (2.24), (2.19), respec- 
tively. 

Recall that for a U(N) gauge group, <Tn is a constant N x N Hermitian matrix that 
commutes with m. We may thus diagonalize 

/ Ai Atv 

where Xi/2n, i — 1, . . . , N, are the real eigenvalues of a , and we order Ai < ■ ■ ■ < A^r- 
A choice of positive roots for G is then 

A,; — A, 



(2.26) 



2vr 



(2.27) 



with i < j. Notice that the Vandermonde determinant then contributes a factor to the 
integrand of (2.25) when rewriting 



dan 



N 



dX 



m n 

i=l i<j\m,i=mj 



Xi — Xj 
27T 



(2.28) 



which precisely cancels the denominator in (2.24). 

We now specialise to the particular gauge theory of interest, namely the ABJM 
theory on iV M2-branes [8]. This is a U(N) x U(N) Chern-Simons gauge theory with 
Chern-Simons levels (k,—k) for the two gauge group factors, two chiral matter fields in 
the bifundamental representation (N, N), and two chiral matter fields in the conjugate 
(N, N) representation. More precisely, this is the low energy worldvolume theory on 
iV M2-branes transverse to C 4 /Zfc, where the acts with weights (1,1,-1,-1) on 
the four complex coordinates. The R-charges/scaling dimensions of the 4 chiral fields 
all take the canonical value of A = |. We may thus introduce eigenvalues Xi, Xi, 
i = 1,...,N, for the two gauge group factors, and correspondingly matrices m, m 
specifying the flat connections for each copy of U(N). Note that the weights for the 
bifundamental representation (N, N) are 



Xi — Xj 
2vr 



(2.29) 
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with minus this for the conjugate representation. Thus the partition function for the 
ABJM theory on L(p, 1) is 



1 f N 

m,m 

j J 2 sinh 
Y[ 2 sinh 

i<j 

n 



2tt" 27T 



cxp 



ik 

Anp 



N 



N 



A 



(27r) a x;w-^ 2 ) 



,i=i 



i=l 



A, — Aj + 27Ti(m i — in j ) 
2^ 

Aj — Aj + 2-n\[jhi — rhj) 



2p 



2 sinh 



2 sinh 



Aj — Aj — l-Kiivrii — ra j) 
2p _ 

Aj — Aj — luxirhi — rhj) 
2p 



x 



x 



'■j 



pmi-rrij 
P 



\i - Xi 



2tt 



(2.30) 



where we have defined 



n 

£=0 



, 3 , • \ 



(2.31) 



The latter is precisely the contribution from one (N, N) chiral field, and one (N, N) 



field, and these correspond to the respective factors in the product (2.31). As before 



the notation N K (£) means the number of half-integers n 6 {— |, — 4 + 1, •••,!} such 
that 2n = k mod p. The square at the very end of ( 2.30[ ) then accounts for the fact 
there are two of each type of bifundamental field. 



2.6 Matter potentials 

The matter potentials 



v;(z) 



(2.32) 



where P£{z) is defined by (2.31), play an important role in the dynamics of the ma- 



trix model (2.30). A general discussion of these potentials, which in general involve 



polygamma functions, may be found in appendix |A| In particular, the products in 



(2.31) are divergent and must be regularized, and we do this using zeta function reg- 



ularization. The resulting functions simplify somewhat in particular cases. In this 
subsection we give a few examples, for low values of p. 

Recall that for fixed p > 1, we have < k < p. The (regularized) matter potentials 
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for small p then simplify to 



V°{z 
V 2 \z 

V 2 \z 

V 3 °(z 
V 3 \z 



— log [2 cosh(7T2;)] 
1 



-- log[2 cosh(zTr)] + 2z cot" 1 e nz + - [Li 2 (-ie~ OT ) - Li 2 (ie" 

2 7T 

-- log(2 cosh(zTr)) - 2z cot" 1 e nz - - \Li 2 {-ie- Z7T ) - Li 2 (ie 

2 7T 



log 



2cosh 2 (vr2/3) 



cosh 71 z 

V^iz) = -log[2cosh(7rz/3)] . 



(2.33) 



The reader is referred to the Appendices [A] and O for more detail. 



3 The large N limit: numerical results 

Our aim in the remainder of the paper is to compute the M-theory limit of the ABJM 



partition function (2.30), which means fixed Chern-Simons level k and N — > oo. The 



partition function (2.30) is an extremely complicated object, and to gain some intuition 



we will begin in this section by solving the matrix model numerically for large values 
of N. We will do this by extending the saddle-point methods of [9] to the present case. 
The behaviour is simple enough to suggest an ansatz for the eigenvalue distribution, 
precisely as for the ABJM model on S 3 studied in [S], which in section [4] we then 
analytically show reproduces the numerics. Moreover, this analytic result agrees with 



the expected large iV gravity dual result for the free energy of Fj 



M— theory 



3p 



N 3 / 2 . 



3.1 General discussion 



The matrix model partition function (2.30) of iV M2-branes on a Lens space L(p, 1) 
has the following form 

N 



m,m 




JJdAidAi 



(3.1) 



i=l 



where recall that m = diag(mi, . . . , m^) has entries mi < m 2 < ... < with 
< rrii < p, and specifies the flat connection for the first U(N) gauge group, while 
tilded quantities refer to the second U (N) gauge group. The basic idea is that when the 
number N of eigenvalues Aj is large, each contribution Z m ^ can be well approximated 
in the saddle-point limit by Z = e~ F , where the free energy F is an extremum of 
-Pm,m(Aj, Aj) with respect to Aj and Aj. 
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Given F(Aj,Aj) (in what follows, we suppress the lower indices m, rn) the saddle- 
point equations are 

dF dF 

Mr ' air ' ^ (32) 

The extremum of the free energy is then given by F(A°,A°), where A°, A° are the 
solutions of the saddle-point equations. This then gives the leading contribution to 
Z m ,m at large N and fixed k. As for p — 1, it will turn out that the saddle-point 
solution has complex eigenvalues, which means we deform the real integrals over Aj, Aj 



in (3.1) into the complex plane. 



Even though highly non-trivial, the equations (3.2) can be solved numerically. It is 
convenient to view these equations as describing the equilibrium configuration of 2N 
particles, whose two-dimensional coordinates are given by the complex numbers Aj and 
\i. This equilibrium configuration can be found by introducing a "time dependence", 
so that Aj, Aj — > Aj(t), \i(t), and writing down equations of motion for Aj(t) and Aj(t) 
such that their solutions approach the equilibrium configuration for late timesj^] 

d\i _dF_ d\ _dF_ 

dt d\ dt ok' 1 ' ' 

In the following we will solve these equations numerically. From this we will extract 
generic behaviour that will lead to a corresponding analytic computation in section |4j 

3.2 Flat connection dependence 

A new ingredient in the partition function on Lens spaces, with respect to the case on 
S 3 , is the sum over different flat connections labelled by m, m. We are interested in 
the large N limit, and in this limit we expect 



e -^ m ,* _^ e -^M- thoory = Z M -theory • (3.4) 



In the supergravity approximation to M-theory, we are computing the log of the par- 
tition function in the large N limit, and we are interested in the leading term only. 



Hence, even though we have the sum of many terms on the left hand side of (3.4), 
we expect only certain terms to contribute. More precisely, we may focus on the con- 
tribution of (m, m) = (mo,m ) with least F mo ^ Q , in the large N limit. Note that 



order for the eigenvalues to go to the correct attractor point as t — ¥ oo, a priori one might 



need to multiply the left hand sides of (3.3) by a complex number. For the case at hand in fact this 
is unnecessary. 
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contributions of e - ^* for other choices of (m, m) do not need to be suppressed with 
respect to that for (mo, mo): if they give a similar contribution, this will simply lead 
to a logarithmic, and hence subleading, correction at large N, since we are taking the 
logarithm to obtain the free energy. 

We have performed a numerical evaluation of -F m ,m, for several values of N,p and 
all choices of m, m. In all cases the leading contribution comes from m = m = 
diag(c, c, c, c) = c ■ lAr X iv, where c is an integer with < c < p. Hence the first 
lesson we draw from the numerical analysis is that we may focus on the specific case 
m = m = diag(c, c, c, c) if we are only interested in the large N limit. This simplifies 
the problem, and its treatment, enormously. Furthermore, note that with these choices 
of flat connection the eigenvalues will respect certain symmetries, discussed below, but 
that this would not be true for any other choice of flat connection. So from now on we 
focus on this caseJ3 

3.3 Numerical plots 

Figure [T] shows the distribution of eigenvalues for the case p = 2 and N = 100. From 
this distribution we can draw several conclusions. First, we see that the eigenvalue 
distribution is invariant under A, — > —A*, Aj — > — Aj. Furthermore, for the equilibrium 
configuration A, and Aj are complex conjugates of each other. To be more precise, we 
find that Aj = — Atv-i+i (with the same for Aj) and A, = Aj. As for the p = 1 case, 
these are symmetries of the equations of motion, so we expect these symmetries for the 
equilibrium distributions as well. As already mentioned, however, these symmetries 
will not be present for more general choices of m, m. 

Other features are that for large values of N the density of eigenvalues is relatively 
uniform, the real part of the eigenvalues grows with N, while the imaginary part stays 
bounded. The numerics are consistent with the real part growing as N 1 ^ 2 , while the 
imaginary part of the eigenvalues stays bounded between — 7r/2 and tc/2. As we increase 
the Chern-Simons level k the slope also increases. The analytic treatment in section [4] 
(after assumptions justified by the numerics) predicts a slope proportional to Vk - in 
fact precisely the same slope as for p = 1. This is also consistent with the numerical 
results - see Figure [2} 

2 Other choices of m, m, presumably important for computing subleading corrections, are briefly 
discussed in Appendix B. 



13 



p=2,N=100 
Im(A) 




Re(A) 



Figure 1: Distribution of eigenvalues for p = 2, N = 100. 
p=l,2,3,N=100 

Im(A) 




Re(A) 



Figure 2: Distribution of eigenvalues for p = 1,2,3 and N = 100. 

4 The large TV limit: analytic results 



As in the previous section, the idea is to compute the partition function (2.30) in a 
saddle-point approximation, focusing on the contribution from m = m = c ■ Inxn, 
which from the numerics we see determines the free energy in the large N limit. As the 
number of eigenvalues N for each gauge group tends to infinity, one has a continuum 
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limit in which one can replace the sums over eigenvalues in the potential by integrals. 
In particular, one can then separate the interactions between eigenvalues into "long 
range forces," for which the interaction between eigenvalues is non-local, plus a local 
interaction. A key point is that these long range forces automatically cancel. An 
appropriate ansatz for Aj will then lead to a simple local action for the eigenvalues, 
which may be solved in the saddle-point approximation exactly in the large N limit. 
This analytic result may then be checked against the numerical results, and we find 
excellent agreement. We will also comment on the relation to the gravity dual. 

4.1 Long range forces 

Let us focus first on the long range forces, which come from the leading terms in an 
asymptotic expansion of the sinh and matter potential Vp =0 (z) in ( |2.30[ )f|ln the former 
case we define 

[log 2 sinhz] asymp := z sign (Re z) . (4.1) 
The point here is that for Re z > we have the series 

XI j 

log [2 sinh z) = z - -e" 2fe . (4.2) 

while for Re z < we have 

log [2 sinh z] = i7c-z-J2 - p e 2lz ■ ( 4 - 3 ) 



We shall see momentarily that the constant in term in (4.3) does not contribute to 



the long range force computation, which is why we omit this constant in the definition 



(4.1). The sums of exponential terms in (4.2), (4.3) will be of relevance momentarily. 



The matter potentials Vp(z) depend in a complicated way on p. The relevant asymp- 
totic expansions are discussed in appendix [Cj In particular, this leads to 



K(z)r mp := sign (Re.) . (4.4) 

We then have the following general form of the expansions for VP(z) 

V°( Z ) = [y p (2)] aSymP + 7r.sign(Re.) ^ Qe -^sign(Re,) + dee -^si g n(Rez) ^ 

p p 
3 Recall here that since m = m is proportional to the identity matrix, we have K = in all cases. 
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for appropriate constants eg, dg depending on p. 

In the large N limit we then take a continuous limit, in which sums over % = 1, . . . , N 
become Riemann integrals 



N 



N 

E 

i=i 



p{x) dx. 



(4.6) 



The numerical results of section [3] then suggest we make the following ansatz for the 
eigenvalues 



X(x) 



N a x + iy(x) 



X(x) = N a x-\y{x) 



(4.7) 



where a > 0, and these formulae are understood to be correct to order N e , for some 
e > 0. Notice we have deformed the real eigenvalues of the Hermitian matrix a into the 



complex plane in (4.7), anticipating a complex saddle point, and that the function p(x) 



describes the density of the eigenvalues. Also recall that \ <->■ Aj is a symmetry of the 



system - in (4.7) we have imposed that the solution is invariant under this symmetry, 



which is again supported by the numerical results. 

The long range forces are then, by definition, determined by the leading asymptotic 



terms in the potential. Substituting (4.1), (4.4) into the logarithim of the partition 



function (2.30) and taking the continuum limit (4.6), we obtain 



- F, 



asymp 



N 2 



p(x) dx 



p(x') dx' sign(x — x) 



2p 



(X(x) - X(x')) 



i , ~ , , w s/s ^ M%) — x(x') 

+ -(AM - A(x)') - 2 x LL_LJ 



(4. 



Here one sees from (2.30) that one substitutes z = (X(x) — X(x'))/2n into (4.4). The 



factor of 2 in the last term of (4.8) accounts for the two copies of chiral fields in the 



(N, N) and (N, N) representations. Also notice that the original sum in the vector 



multiplet contribution to (2.30) is over i < j, which means x — x' < in the continuum 



limit. In writing (4.8) we have simply extended this to a sum over % > j by replacing 



Xi — Xj by Xj — Aj. It is then straightforward to see from the ansatz (4.7) that all the 



terms in (4.8) cancel: the real parts simply cancel inside the square bracket, while the 



imaginary parts contribute zero on using the anti-symmetry under x •<->■ x' implied by 
the sign(x — x') term. Thus -F asymp = 0, and the long range forces indeed cancel for 
L(p,l). 
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4.2 Local action 



It follows that only the exponential sums in (4.2), (4.3), (4.5) contribute to the partition 



function in the large N limit. Since the latter depend on sign (Rez), which is equal to 
sign(x — x'), we first split the double integrals as 



N 2 



p(x) dx 



p(x') dx' 



N 2 



+N 2 



p(x) dx / p(x') dx' 

lax /"^max 

p(x) dx / p(x') dx' , (4.9) 



so that x — x' > for the first term on the right hand side, while x — x' < for the 
second term. We will then apply the general formula 



dx'e-? Na ( x -^ f(x,x') 



/3N C 



-0N a (x-x') 



f(x,x') 



(3W 



„n d 



dx , e -PN°>( x - x >) J^f( x J) (4.io) 
dx' 



which follows trivially from an integration by parts. The first term on the right hand 
side is simply j^f(x,x), plus a term which is exponentially suppressed in the large 



N limit. The formula (4.10), with a similar formula applying for x — x' < 0, amount 
to the representation 



c 



6(x) = lim -e" c|x| , (4.11) 

c— >oo 2 

thus reducing the integral over x, x' to an integral over x, in the large N limit. Applying 



this to the sums over exponentials in (4.2), (4.3), (4.5) is a straightforward task. For 



the vector multiplet and matter multiplet contributions, we obtain to leading order 

4pir 2 



^vector 
-F mattcr 



6 



-N- 



2-a 



p{x) 2 dx + o(N 



2-a\ 



(4.12) 



d 



2-a I i \2 , I °£ 



p(x) 2 dx<{ ^2 -jy( x ) sin [fy( X )] + p COS 



x 



(4.13) 



The term in curly brackets is denoted J p [y(x)} in Appendix [Cj and may be evaluated 
by Fourier summation to give 

OX y{x) 2- 



- F, 



matter 



IN 



2-a 



p(x) dx 



24\ P ~p) + 2p 



+ o(N 2 ~ a ) .(4.14) 
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Combining with (4.12), we thus obtain the leading order result 

jl-a 



one— loop -^vector 



]\T2—a rxm^x 

+ ^matter = / p{xf ~ Ay{ X f\ dx + o(iV 2 -«) (.4.15) 

P J Itnin 



xy(x)p(x) dx + o(N 1+a ) . (4.16) 



It remains the add the contribution of the classical terms in (2.30). This is a trivial 
modification of the p = 1 computation, the only difference being the factor of 1/p = 
Vol(L(p, l))/Vol(S 3 ) : 

kN l+c 

^classical = 

pn 

The total free energy action is then to leading order 
pn 

As for the case of the round sphere with p — 1, non-trivial saddle points will require 
both terms to be of the same order, so that a = | and hence A(x) = N l / 2 x + iy(x). 



/ xy(x)p(x) dx H / p(a;) 2 [ff 2 - Ay(xf] dx .(4.17) 

" ^min ^ *^ ^min 



Remarkably, we see that the action F in (4.17) is simply 1/p times the action for 
p = 1 in reference [9]. In particular, the saddle point equations derived from (4.17) 
are identical to those in reference [H], which allows us to simply write down that the 
density p(x) is constant 

p(x) = — -=- , (4.18) 
and the imaginary part of the eigenvalues y(x) is linear 

with — x min = x min = TTy/2/k, so that y{x) £ [— Of course, this is perfectly 
consistent with the numerical results in Figures [T] and [2] The dependence on p only 
enters in the free energy F evaluated on this saddle point solution, which is 



N^-^- + o(iV 3 / 2 ) . (4.20) 
3p 



Again, this is consistent with the numerics. 



The formula (4.20 ) is expected from the supergravity dual solution AdS^/Zp x S /Z^, 
since the quotient by Z p simply divides the overall supergravity action by p. The only 
slight subtlety here is that AdS 4 /Z p has a Z p orbifold singularity at the "centre". In 
principle there might exist degrees of freedom at this singularity which then contribute 
to the leading order large N free energy, but the field theory result we have obtained 
implies this is not the case. 
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5 Outlook 



In this paper we considered the large N limit of the partition function of N M2-branes 
on the Lens space L(p, 1). Some open problems include: 



The partition function (2.30) is valid for all N,k,p, and it would be interesting 
to study this more generally, for example at finite N, or in the 't Hooft limit in 
which N/k is held fixed. 

One might also consider squashed Lens spaces, for which there are supergravity 
dual solutions 



• Another interesting open question is whether these theories have a description in 
terms of a Fermi gas, as for the ABJM theory on S 3 [13J. This may be a useful 
method for computing subleading corrections. 

The generalization of these results to more general Lens space L(p, q) will be addressed 
in a forthcoming publication. 

Acknowledgments 

L. F. A and J. F. S. would like to thank the Isaac Newton Institute for hospitality 
during the completion of this work. J. F. S. is supported by a Royal Society Research 
Fellowship. 

A Computation of potentials 

In this appendix we present analytical expressions for the infinite products that enter 
into the partition function for the Lens spaces L(p, 1). These infinite products have 
the following form: 

where k = 0, l,...,p — 1 and N K {£) denotes the number of integers m = {—£,—£ + 
2, ...,£ — 2,£} such that m = k mod p. For computing the free energy the log of this 
product is relevant. Hence we introduce the potentials V*\ 

V; := logP p K . (A.2) 
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Let us explain in detail how to obtain the potential for p = 1. This result is already 
known [2], but it is instructive to recover it. In this case we have k = and the 
potential reduces to 

KM = f>+W^^-^4^) • ( A - 3 ) 

This sum is divergent, hence in order to compute it we need to regularize it. A standard 
procedure is to take the derivative of the potential and perform the sum. We obtain 

V[(z) = -Trtanh(Trz) . (A.4) 

Of course, in taking the derivative we are dropping an additive constant (which could 
be infinite). Integrating back we find 

V\(z) = — log cosh (tcz) + c . (A. 5) 

The integration constant can be fixed by zeta function regularization. This is explained 
in detail for instance in appendix A of [14J. The constant c is defined as the value of 
Vi(z) at z = 0. We obtain 

2 C = E^ +1 ) lo g(^f) • ( A - 6 ) 
We will compute this divergent sum by using (^—function regularization. Let us define 

, / X ^ ( £ + 1 1 + M /A^ 

Cz(s) = S w " w ■ ( ' 

Hence the quantity we wish to compute is just —(' z (0). These sums are by definition 
zeta functions, and their generalization, Hurwitz zeta functions 

Sums with factors of I in the numerator are easily obtained, since a factor of £ + a in 
the numerator can be absorbed by a shift s — > s — 1. For the particular case at hand 
we obtain 

Us) = -(2 s -l)C(s). (A.9) 
Hence — Cz(0) = — ^p- This implies c = — log 2, and the following result for V\\ 

V^z) = -log(2cosh(7T2)) , (A.10) 
which coincides with the known result. 
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A. 2 General p and k 

Let us start by introducing some more notation: 

m == log (£f±£) . (A.11) 

The potentials will be given by the sum of two terms 

V p K (z) = U;(z) + U^ K (-z) . (A.12) 

By working out some explicit examples, one can convince oneself that the general form 
of each contribution U is as follows 

oo 

Ufa) = + s if(P* + 1) + - + s p -xf{p£ + (p - 1))] , (A.13) 

1=0 

where so, ...,s^_i depend on p,n and £ The important fact is that they are always 
of the form = etj + can be computed in two steps. First we use the 

intermediate result 

°° 1 r /c + iz\ 
^(a + 6£) log(p£ + c + iz) = - (bc-ap + izb) log T 

P l \ P J 

- bp ^(-v(^±h)] , (A.14) 



p 

where tp(~ 2 ^ is the polygamma function, and we have dropped a term that can later 
be fixed (for the final result) by using zeta function regularization. All our expressions 
are the sum of such building blocks. In order to assemble the correct building blocks, 
we just have to compute Sj = aj + bit for the fixed value of p, k that we are interested 
in. Finally, once we have computed the Sj, we can compute the correct integration 
constants by using zeta function regularization, as shown above. It is straightforward 
to write a Mathematica code that computes the final potential Vp{z) for any choice of 
p, The point is that since is at most linear in £, we can compute by looking at 
the terms with < I < 2p — 1 . 

A. 3 Some explicit examples 

Below we present some explicit results that are used in the numerics. Even though the 
general answer depends on polygamma functions, for some cases the final expression 

4 Thc code is available upon request. 
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can be simplified: 



— log [2 COsll(7Tz)] , 

-- log[2 cosIi(ztt)] + 2z cot^ 1 e wz + - [Li 2 



- Liofie" 



--log(2cosh(z7r)) -2zcot~ 1 e KZ - ^ [Li 2 (-ie" 27r ) - Li 2 ( 
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7T 



log 



2 cosh 2 (7^/3) 



V?{z 

V?(z) -- lot; 2 cosh ( - — } | ■■■■ 2.: oo1 1 c"- + - |Lio(-io 

V 3 \z) = V-l^z) = -log[2cosh(7rz/3)] . (A. 15) 

A general formula is given in appendix C for the case k — 0. Some comments are in 
order. We see that the sum of potentials over all k for fixed p satisfies a completeness 
condition 



p-i 



(A.16) 



K=0 



This is of course expected, since fixing k projects over certain terms in the sum giving 
V\. Another comment is that from the structure of the sums we expect Vp(z) to be an 
even function and Vp(z) = V p p ~ K (— z). 



B A wave of eigenvalues 

In the body of the paper we have shown that in the large N limit we can focus on the 
case m = rh = c ■ Inxn- Furthermore, we have analyzed numerically the distribution 
of eigenvalues for this case. One can use the numerics to analyze the eigenvalue distri- 
bution for other choices of m and m. These will presumably be important if one wants 
to compute sub leading corrections to our result. An interesting distribution of eigen- 
values is obtained ifm = m^c-lAr X Ar - see Figure [3} We have shown the eigenvalue 
distribution for p = 2 and N = 100, for the cases m = m = (0, . . . , 0, 1, . . . , 1), with 
n zeros and N — n = 100 — n ones. From left to right, top to bottom, we show the 
eigenvalue distribution for n = 3, 10, 20, 30, 50, 70, 80, 90 and 97. The distribution of 
eigenvalues is reminiscent of a wave moving from left to right, with the location of the 
kink at the boundary between the group of zeros and the group of ones. 

Equivalently, we see that for p = 2 and m = m the eigenvalues distribute in two 
segments. The length of the segments is equal to the quantity of zeros and ones, 
respectively. The numerics seem to suggest that when the number of zeros and ones 
is "macroscopic" (i.e. of the same order as N) there is a finite "jump" between the 
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Figure 3: Distribution of eigenvalues for N = 100, p = 2 and n — 3, 10, 20, 30, 50, 70, 
80, 90, 97 from left to right and top to bottom. 

segments. This feature is also present for other cases. For instance, Figure [4] shows the 
eigenvalue distribution for p = 3, N = 100 and m = rh = (0, . . . , 0, 1, . . . , 1, 2, . . . , 2), 
with 30 zeros, 30 ones and 40 twos. It would be interesting to understand whether the 
finite "jump" between segments is really there for large N or an artifact of iV being 
not large enough. 

C Asymptotic expansions 

In this appendix we present the asymptotic expansions for the potentials found above. 
We will present a detailed analysis for k — 0, since, as discussed above, this is enough 
to compute the free energy in the large N limit. 

Proceeding as explained in Appendix A we find the following expressions for p odd 
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p=3, m=(0„ 



,0,1,-,1,2...2) 
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40 







-0.5 - 



-1.0 - 



Figure 4: Distribution of eigenvalues for N = 100, p = 3 and m = m = 
(0, . . . , 0, 1, . . . , 1, 2, . . . , 2). We see that the eigenvalues distribute in three segments. 



and even respectively: 



9z Vp =0 dd 



d V 

w z v p=even 



• P-1 
1 

V 



n=0 



1 + 2n + 2iz 
2p 



■0 



n=0 



(2z + i(2n + 1 -p))V> 



1 + 2n - 2\z 
2p 

1 + 2n + 2i^r 



2p 



+(2«-i(2n+l-p))V' 



1 + 2n - 2iz 
2p 



(C.l) 



For p odd these expressions can be given in terms of trigonometric functions, but the 
present form is more uniform. 

Now we would like to compute the asymptotic expansions of such expressions, when 
the real part of z is very large. We have the following for \z\ — > oo and arg(z) very 
close to 7r: 

1 °° B 1 
tftz) = log(z) - - - J2 ^ + >0 ^t(vrz) - 1) , (C.2) 



2z ^ 2kz 2k 2 v v ' ' ' 

k=l 

where B 2 k are the Bernoulli numbers. The coefficient in front of i7r(i cot(7rz) — 1) is 
actually one if | arg(z)| > ir, and is zero otherwise. For the case of real functions, as 
we are considering, these two average to 1/2. 

It is now easy to compute the asymptotic expansion, substituting the expansion for 
ijj into the expression for d z V p . For each case we obtain: 
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C.l p odd 

d z V p=odd = ~2---J2^i^- (C3) 
This gives the following expansion for lp =c ,dd: 



oo -2l 2 < 



V p=odd = -2-z + V 6 % e , ■ (C4) 
P ^cos(f) 

Quite remarkably, zeta function regularization implies a value for V(;z = 0) such that 
the asymptotic expansion doesn't have a constant term. 
Given an expansion of the form 



V = zir (co + + X> 



' e e-* zi , (C.5) 
we have seen in the body of the paper that the contribution relevant at large N is 

J = E ( j y sin (^) + % <x»(*y) + j cos ( £ y)) ■ ( c - 6 ) 



e 

For the present case we have 



p ~ cos(fy) 2 p 2 -3 2 y 2 



^=odd - 2 Z. £2cos(£ p - 24p + 2p' 1 J 

where we have assumed that y lies in the range [— ir/2, 7r/2]. This is justified by the 
numerics. 

C.2 j9 even 

For p even it is convenient to focus on the contribution for a fixed value of n in the 
general expression for d z V P:1 . We obtain 

\n ri - T ,^7r(2n+l) ^ 



-e p 



P" P ^ P 



+ 2(-ir(l + 2n-p), Eco / ( 2n+ l), A _ 2; „ (C8) 



V "VP 

The contribution from this term to J, which we denote J\ n , can be computed as above. 
We obtain 

/ (1 + 2n - + 4n 2 - 4n(p - 1) + 2p(p - 1)) 2 p - 1 - 2w 2 

J| " = 24? * + 2 P 2 y 

(C.9) 
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Summing over n from zero to p — 1 and using that p is even we obtain 

J P =eve„ = -^rr + - . (CIO) 

This has exactly the same form as for p odd! We thus arrive at the following result, 
valid for all values of p: 

2p 2 -3 2 , ?/ 2 



J v = — —71 + ^- . (C.ll) 

p 2Ap 2p y ' 
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